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Abstract 

We present in a unified and self-contained manner the coordinate-free approach to spherical 
harmonics initiated in the mid 19th century by James Clerk Maxwell, William Thomson and 
Peter Guthrie Tait. We stress the pedagogical advantages of this approach which leads in a 
natural way to many physically relevant results that students find often difficult to work out using 
spherical coordinates and associated Legendre functions. It is shown how most physically relevant 
results of the theory of spherical harmonics - such as recursion relations, Legendre's addition 
theorem, surface harmonics expansions, the method of images, multipolar charge distributions, 
partial wave expansions, Hobson's integral theorem, rotation matrix and Gaunt's integrals - can 
be efficiently derived in a coordinate free fashion from a few basic elements of the theory of solid 
and surface harmonics discussed in the paper. 
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I. INTRODUCTION 



Spherical harmonics present themselves in many branches of physics and applied math- 
ematics such as electromagnetism, quantum mechanics and gravitation theory. In most 
modern textbooks dealing with topics of mathematical physics^, the theory of spherical 
harmonics is built upon their representation in spherical coordinates. As is well known, 
in this approach one starts by expressing Laplace's equation in spherical polar coordinates 
(r,9,ip), where r is the distance to the origin and 9 and ip are the polar and azimuthal an- 
gles, respectively. One then tries a solution in separated variables that yields the spherical 
harmonics in the form Yi m (9, ip) = Ci m e im,fi P™ (cos 9) (I = 0, 1, —l<m< I), where C/ m 's 
are normalization constants and the P z m (/z)'s are the associated Legendre functions defined 
by 2 l l\P l m (fj,) = (-l) m (l - fi 2 ) m / 2 d l+m (ft 2 - l) l /dfi l+m . However, this representation often 
leads to cumbersome derivations of the main results needed in physical applications, so that, 
from a physicist's point of view, it may be desirable to have a more compact approach in 
which mathematical manipulations become more direct and, perhaps, also more guided by 
physical insight. This is also the case with other special functions of great physical relevance 
such as Bessel's functions or ellipsoidal harmonics. 

The first representations of spherical harmonics freed from spherical coordinates were 
those given by Thomson and Tait 3 in their Treatise on Natural Philosophy and by Maxwell 4 
in his Treatise on Electricity and Magnetism. The principle of their method is that, having 
given any solution of Laplace's equation, other solutions may be obtained by differentiating 
the given solution any number of times in the directions of the cartesian axes x,y and z or, 
as Maxwell proposed, in any directions whatsoever. In particular, by applying this method 
to the simple solution 1/r, with r = x 2 + y 2 + z 2 , all the system of spherical harmonics 
may be generated^. Maxwell 4 - also contributed with an important integral theorem which 
apparently has not deserved by itself much attention in the literature. In fact, after it first 
appeared in Maxwell's treatise, it has almost invariably been presented as a corollary of 
a more general - but also more complicated - integral theorem first introduced by Niven- 
and generalized afterwards by Hobson^. Nevertheless, we consider Maxwell's theorem as an 
essential ingredient in our approach and, therefore, have preferred to follow here a different 
route that leads to Hobson's theorem via Maxwell's theorem. 

The aim of this paper is to show how the contributions due to Maxwell and to Thomson 
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and Tait lead to a unified, self-consistent and spherical coordinate-free approach through 
which most of the properties of spherical harmonics needed in usual physical applications 
can easily and naturally be derived. We hope that this paper helps students and researchers 
to understand and work out by themselves many of the results they need when dealing with 
specific problems involving spherical harmonics, thus avoiding much of the trouble of con- 
stantly having to look them up in the appendices of the relevant literature or in specialized 
mathematical handbooks^. For these purposes, we have based our presentation on the sim- 
ple properties of three basic elements of the theory of solid and surface harmonics: a) the 
elementary solid harmonics of the form (b • r) l -l being a nonnegative integer and b a null 
vector-, b) Maxwell's harmonics, and c) Maxwell's integral theorem. We have introduced 
these basic concepts in Section II and have applied them in the rest of the paper to pro- 
vide systematic, coordinate-free derivations of some physically relevant results such as the 
construction of the standard set of spherical harmonics Yi m (Section III), recursion relations 
(Section IV), Legendre's addition theorem (Section V), surface harmonics expansions (Sec- 
tion VI) - with applications to the method of images and multipolar charge distributions-, 
partial wave expansions and Hobson's integral theorem (Section VII), rotation matrix (Sec- 
tion VIII) and Gaunt's integrals (Section IX). Finally conclusions are presented in Section 
X. 



II. GENERAL PROPERTIES OF SOLID AND SURFACE HARMONICS. 

We define a solid harmonic of the type Hi or, more briefly, a Hi-harmonic, as any solution 
of Laplace's equation, 

V 2 = O, (1) 
which is a homogenous polynomial of degree I in the (cartesian) variables (x, y, z): 

H l (r)= CWV* 7 , (2) 

where r = xe x + ye y + ze z is the position vector and constants C a/ g 7 are generally complex. 
As can be checked by direct substitution, to every inharmonic there corresponds a second 
solution of ([I]) of the form 
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which we will define as a solid harmonic of the type VJ or, more briefly, Vi-harmonic. Note 
the following behaviors of Hi and V\ for small and large values of r 

r^O: #,(r)->0, VJ(r) oo, (4) 

r -»■ oo : iZ",(r) oo, VJ(r) 0. (5) 

According to (jlj), if /-harmonics and V^-harmonics are also known as regular and irregular 
solid harmonics 9 -, respectively. 

A surface harmonic of degree /, YJ, is defined as the function obtained from a inharmonic 
or, equivalently, a VJ-harmonic according to the rules 

tf(r) = #,(r)/r' = r m V^(r). (6) 

Since Hi is a homogeneous polynomial of degree I in the variables x, y and z, it is clear that 
the value of Yi at any point of space, r, is the same as that at the point of the unit sphere 
e r = r/r: Yj(r) = Yi(e r ). Also, note that the values of Yj, if; and V\ coincide on the unit 
sphere. 

It can easily be shown^^ that the number of independent Hi - harmonics is 21 + 1 which, 
of course, must also be the number of independent functions Vi and Y\ by fl3]) and fl6]). In 
effect, to a given value of the exponent a (0 < a < I) in a homogeneous polynomial of degree 
I such as (jSJ) there correspond I — a + 1 possible values of the exponent /3, since the value of 
7 is fixed by 7 = / — a — (3. Thus, the maximum number of monomials of the form x a y^z 1 
in (j2J) will be ^^ =0 (^ — a + 1) = (I + + 2)/2. Since V 2 i?z is a homogeneous polynomial 
of degree 1 — 2, it will therefore contain a maximum of (/ — 1)1/2 terms whence the condition 
V 2 Hi = is equivalent 1(1 — l)/2 equations between the constant coefficients C Q( g 7 in (J2J). 
The number of independent constant remaining is thus (I + 1) (7 + 2)/2 — 1(1 — l)/2, or 2/ + 1 . 
Therefore, any solid or surface harmonic of degree / can be expressed as a linear combination 
of 21 + 1 independent harmonics of the same class. 

Observe that, unlike functions Hi or Vi, surface harmonics do not satisfy Laplace's equa- 
tion, but instead 

V 2 Yi + fcll^ = 0, (7) 

as is readily verified by direct substitution of ([6]) into ([I]). From (J7]) it is easy to obtain the 
orthogonality property^ for two given surface harmonics, Yi and X n . In effect, if we substract 
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the equation for Y\ multiplied by X n from that for X n multiplied by Yi and integrate over 
the volume of a sphere of radius a we obtain 

[1(1 + i) _ n (n + i)] I dV^YtXn = - f dS (X n e r • VY~ ; - Y lGr ■ WX n ) = 0, (8) 

J r<a T J r=a 

where the surface integral has been obtained through Gauss's theorem on using the identity 
X n V 2 Yi — YiV 2 X n = V ■ (XnVYi — YiVX n ). Since surface harmonics do not depend on r, 
the surface integral vanishes and the volume integral in the LHS of (jHJ) can be transformed 
into one on the unit sphere by writing dV = r 2 drdQ, dQ being the element of solid angle. 
Carrying out the trivial integration over the radial coordinate © yields 



/ 



dQYiX n = if l^n. (9) 



An interesting relation between the solid and surface harmonics satisfying ()6]) is obtained 
if one considers the function 
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- [a- l+1 H[(r) - a* +2 V;(r)] , (10) 



21 + 

where a is any positive constant. On using ()6]) to express Vi and Hi in terms Yj and taking 
into account that e r ■ VYj = it is seen that ^ vanishes at r = a and satisfies the identity 

<9# 



1 dr 



[er-V*] r=a . (11) 



The physical meaning^ of this identity is that a~ l+1 Hi(r)/(2l + 1) and a l+2 Vi(r) / (21 + 1) 
yield the potential for r < a and r > a, respectively, due to an arbitrarily thin sphere 
of radius a with a surface charge density distribution given by Yi(e r )/(47r) - the jump at 
the surface in the normal derivatives of the potential being represented by ffTTT) -. This 
expression is useful when one considers the integral over the sphere r = a of the product of 
Y\ times any regular function $ that satisfies V 2 $ = inside the sphere. Using the identity 
V • ($W) = V ■ (*V$) + $V 2 ^ - ^V 2 $ with V 2 $ = 0, Gauss's theorem yields 

/ dS^Yt = I dV<$>V 2 ^ = r / dV$V 2 V h (12) 

J r=a Jr<a ^'~\~*-Jr<a 

where we have taken into account in ffTUl) that Hi is regular for r < a and V 2 Hi = 0. In f[T2"j) 
\7 2 Vi must be considered as a generalized function due to its singularity at the origin. 
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A. The simplest solid harmonics 



The simplest if /-harmonics are polynomials of the form (b • r) , where b is any constant 
(complex) null vector, i.e., b • b = 0, so that 

V 2 (b • r) 1 = IV ■ [b (b • it)' -1 ] = 1(1 - l)(b • r)'- 2 b • b = 0. (13) 

It can be shown that any if;-harmonic can always be written as a linear combination of 
the form 

ff,(r)=A 1 (b 1 -r) l + A 2 (b 2 -r) J + ..., (14) 

for suitable constants Ai,A%,... and null vectors bi,b2, ... whence any property that can 
be easily found out for the simplest solid harmonics can be extended by linearity to any 
general To prove (fl4|) . let us consider the particular set of null vectors of the form 
b = i cos m e x + i sin-ue y + e z (0 < u < 2n), or (b • r)' = (ixcosu + iysinu + z) 1 . On 
expressing the trigonometric functions in polar form and binomially expanding we find for 
any given value u s in (0, 2n) that 

i 

(ixcosu s + iy smu s + z) 1 = ^ hi {k '\x, y, z)e ik ' Us , (15) 

k'=-l 

where the coefficients of the trigonometric series, i.e the homogeneous polynomials 
h^ k \x,y,z), must constitute a set 21 + 1 independent solid harmonics. By taking the 
set of points u s = 2irs/(2l + 1) (s = 0,...2Z) and using the orthogonality property 
Eio e- 2nis( - k - k 'V( 2l+ V = (21 + l)8 kkl we can easily invert to obtain 

h^(x,y,z) = j^e 2 ^'^ Uxcos \^\ +zysin \^\ + z\ , (16) 

whence each of the 21 + 1 independent solid harmonics h^ k \x,y,z) can be expressed as a 
linear combination of solid harmonics of the form (bj ■ r) 1 . Since every if^-harmonic is a 
linear combination of the /i/^'s, equation ( fl4l) immediately follows from ( fTBT) . 

An important application of fTl4|) results from the analysis of the repeated action of the 
differential operator b • V - b being a constant null vector - on any spherically symmetric 
function F(r). In effect, since 



and 

b • V(b ■ r) n = n(b ■ r) n_1 b ■ V(b ■ r) = n(b ■ r) n_1 b ■ b = (18) 
for any integer n > 0, on applying / — 1 times (I > 1) the operator b ■ V to (fT7|) we obtain 

(b.V)<F = (b-r)<(l£)V. (19) 

Equation ffl"4]) then implies 

fli(V)F(r) = fl,(r) F(r), (20) 

where Hi(V) is the differential operator obtained by substituting each cartesian component 
of r in Hi(r) by the corresponding one of V. 

As an illustration of (l20l . let us consider the integral (to be used in Section VI on partial 
wave expansions) 

J dn,#,(q)e q - r = Hi(V) J dtt q e* T , (21) 

where q = q x e x + q y e y + q z e z is a vector of the form and q = qe q - q being a (possibly 
complex) isotropic scalar and e q a real unit vector -, and dfl q = sin 9 q d9 q d(p q is the element 
of solid angle around the point e q of the unit sphere. The equality in (121]) easily follows 
from the properties of the derivatives of the exponential function taking into account that 
Hi(q) is a homogeneous polynomial. Introducing fi q = cos8 q = e q ■ e r , the last integral in 
(j2IJl yields 

Jdn q ^^J\^^ ^\ 



(22) 



so that, on using ([2"0l . 

/ rffi g ^(q)e q - r = 47r^(r) 

Defining the function 



(--) 


sinh(gr) 


\r dr J 


qr 



(23) 



we can write ( )23|) more briefly as 

y d^^(q)e q - r = 47r^(r)g 2i S/ (^)- (25) 

Note that the series expansion for sj(£) can be easily obtained by applying the operator 
[(l/C)rf/rfC]' = 2' [d/^C 2 ]' to the power expansion of sinhC/C. The fact that this series only 



contains even powers of ( makes possible to compute (121]) for the case in which the vector r 
is replaced by the gradient operator V. In effect, it is easy to convince oneself that then the 
result must be formally the same as (12"3"1) but with r = y/r ■ r replaced by a/V • V = VV 2 : 

J rffi g ^(q)e q V = 4nH l (r)q 2l s l (qy/v^j , (26) 

where the operator on the RHS must be computed by replacing ( 2 by (qVV 2 ) 2 = g 2 V 2 in 
the series ff24|) which, therefore, yields only integral powers of the operator V 2 . 



B. Maxwell's harmonics and Maxwell's integral theorem 

Maxwell's harmonics'^ are formed by taking I directional derivatives of 1/r in any given 
directions defined by (generally complex) vectors ei,...ej. Thus the function 

m r ) = (e 1 .V)...(e I .V)(l/ r) = |^i, (27) 

where dhi (i = 1, ...I) denotes the infinitesimal (generally complex) displacement associated 
to e i; must be a VJ-harmonic, since it can be seen by inspection that r 2l+1 Vi is a homoge- 
neous polynomial of degree I which satisfies Laplace's equation (since so does Vi). We will 
illustrate this for the first few values of /, and then the interested reader can easily proceed 
by induction^ to prove it for any I. For example, for I = 1 and 1 = 2 we have 

ei"V(l/r) = -^ (28) 

and 

e,.V[e 2 .V(l/ r )] = - 3ei - re --; 5 - r2e '^ (29) 

where it can be readily checked that the homogeneous polynomials in the numerators satisfy 
Laplace's equation. 

Therefore, the function 

is an if^-harmonic, and the function 

Yi(r) 

is a surface harmonic of degree I. Solid harmonics which, except for a constant factor, 
are of the form (j27) and fl30l) will be called Maxwell's solid harmonics of the V\ and Hi 
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,2Z+1 



d '(1/r) 
dhi...dhi 



J+i 



d\l/r) 
dh\...dhi 



(30) 



(31) 



types, respectively, and surface harmonics of the form (13T1) will be called Maxwell's surface 
harmonics. The vectors e« are called the poles of the Maxwell's harmonics, and in the 
particular case in which they are both unit and real, they define I points on the unit sphere 
known as Maxwell's poles^-~—. 

A Maxwell's surface harmonic in which all the ej's coincide with a given vector, say ei, 
is called a zonal harmonic: 



Z l (r)=r^(e 1 .V) l (l/r)=r^^l (32) 

ohi 



Zonal harmonics arise most naturally when one Taylor expands around r = the potential 
due to a unit source localized at r'. For instance, the expansion valid for r > r' can be 
written as 



1 ^ r n (-l)V 



£iA4^(<-V)'(l/r), (33) 



1 1 Z=0 



where r' = r'e' r . Note that all factors r' l /r l+l in the series ( 133]) are multiplied by a zonal 
harmonic of the form (1321) with pole at e^,. On the other hand, the well known- expansion 
in terms of Legendre's polynomials, Pi, is 



00 ji 



1 1 1=0 

whence identifying terms in ( 133|) -( 134|) we obtain: 

(_l)i r l+i 

Pi« ■ e,) = ^ « ■ V)' (1/r) , (35) 

which establishes the relation between Legendre polynomials and zonal harmonics. 

Maxwell's integral theorem results when we substitute Y\ in (JT2l) by a surface harmonic 
of the form (Ell). Then V x = r _ (' +1 )YJ is given by ( 1271) and, since V 2 (l/r) = — 47r<5(r) - where 
5(r) is Dirac's delta function - we have 

Substituting ([36]) into the last integral of ( 1T2]) and integrating / times by parts - taking into 
account that S(r) and all its partial derivatives vanish at the surface r = a -, we obtain 

47r(-l)V +2 d l § 



l=a dS ® Yl 21 + 1 dhx.-.dhi 



(37) 

r=0 

which is Maxwell's integral theorem. The direct proof of ( l37j) given here is the only one we 
are aware of besides that given by Maxwell himself in his treatise. It seems to us than ours is 



somewhat more transparent than his - which is based on the properties singular multipolar 
charge distributions -. 

As an illustration, let us apply ( 13 7p for the case in which $ is any solid harmonic Hi and 
the Maxwell's harmonic is the zonal harmonic Z\ given by (I3"2"j) : 

[_ a dSH lZl = M ~^ +2 [(ex ■ V)^] r=Q . (38) 
In order to compute the RHS of ( 1381) we binomially expand the operator (ei ■ V)', 



(ei-V)-/. 2^ a'\p'\ 1 '\dx<*' dyP'dz^" [ ) 

a'+/3'+7'=« ' y 

and take into account the expression (T5]) for Hi(r). On computing the result at r = we 
obtain 

[(ei-V)^(i-)] r=0 = /J £ C^e^e^eL = (40) 

Q+/3+7=i 

which is a direct consequence of the homogeneity of Hi - and therefore it remains valid if ei 
and Hi(r) are replaced by any constant (generally complex) vector and any homogeneous 
polynomial of degree I, respectively-. Thus 

4:7T(-l) l Va l+2 

dSH ^= 2l + 1 H M- ( 41 ) 



III. THE STANDARD SET OF SPHERICAL HARMONICS Y, 



lin m 



The standard spherical harmonics^- of degree Z, Yi m (r) (m = —I,. ..I), are 21 + 1 inde- 
pendent, normalized surface harmonics such that the m-th one depends on the azimuthal 
coordinate as e imip . The normalization condition is defined as 

d^lYi* m ,Y lm = <W (42) 

As already pointed out, these functions are usually given in terms of the associated Legendre 
functions, P / m (cos6 l ), which can be introduced by separating variables in Laplace's equation 
expressed in spherical coordinates (see also Appendix A). However, we next show how they 
can be generated in an almost trivial fashion by the procedure due to Thomson and Tait£, 
who constructed Maxwell's harmonics by taking directional derivatives of 1/r with respect 
to z and with respect to the complex coordinates 

£ = x + iy = r cos 6e %v and r\ = x — iy = r cos 9e~ %tp . (43) 
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In effect, let us assume for the moment that m > and consider the VJ-harmonic 

d l (l/r) d l ~ m 1 

VUr) = = (-1/2)- (2m ~ 1)H ^ r2m+1 , (44) 

where we have carried out mt/i-partial derivative with respect to rj of r _1 = (?•/£ + z 2 ) -1 / 2 
using the notation (2m — 1)!! = (2m — 1) x (2m — 3) x ...1 with the convention (—1)!! = 1. 
Observe that the presence of the factor £ m automatically yields an azimuthal dependence 
of the form e tmip in (jUJ) . Therefore, the spherical harmonics Y\ m can be obtained from the 
Vim's by just writing 

Y lm = C lm r l+1 V lm , (45) 

where the constants C\ m must be determined from the normalization condition. We can 
readily compute the normalization integral ( 14"2"|) via Maxwell's integral theorem by putting 
$ = r l Y^ Yi = Y im and a = 1 in (13"7j) . which yields 

47r(-l)'|a m | 2 ^(r 2 ' +1 ^ 



= 1. (46) 

r=0 



2/ + 1 dz l ~ m dri 

Note that, since ?" 2Z+1 V^ is a homogeneous polynomial of degree I in the variables £, 77 and 
only the coefficient of the term containing ^ m z l ~ m contributes to the value at r = of 
the derivative in (1461) . Multiplying (144]) by r 2 ' +1 and taking complex conjugates we obtain 

al—m 1 

r 2l+1 VL = (-1/2)™ (2m - 1)!! ^ r^-^ — , (47) 

so that in order to find the coefficient of rj m z l ~ m in (j4T|) we just have to compute the value 
at 77 = and £ = of the factor to the right of r] m , namely 



z 2m 



1 _/ -|V-m (j + m ) ! J-m (A o\ 



Substituting (1481 into ( l4"Tj) we obtain 

r 2* +V * _ (-1)' (2m -!)!!(/ + m)! 

' m ~ 2™(2m)! V + ""' 1 J 

where the dots denote terms which do not contribute to the LHS of ( 1461) . whence 

47r(l + m)\(l-m)\\Ci. 12 



1 1>> 1 



(50) 



(2/ + l)2 2m 

The normalization constant C[ m is determined by (150]) except for a phase factor which we 
will choose as (— l) l+m . This choice is motivated so that the resulting spherical harmonics 
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satisfy the relation (1571) below, as do those used in standard references' 1 ^. Solving (1301) for 
Cim and substituting into (}4"5l) we finally obtain 



*W r ) = \l — m . s T ^- m ^ m - ( 51 ) 



'2Z + 1 (_i)'+^2 m r m d l {l/r) 
4tt ^(Z + m)!(/-^JT^- m ^" 

It is shown in Appendix A that ( 130 reproduces the standard expression— for the IWs when 
expressed in spherical coordinates. 

Until now we have assumed < m < I and found a set of I + 1 spherical harmonics 
given by (130 for m = 0, 1, .../. These functions, together with their corresponding complex 
conjugates Y* m {r) form a complete set of 21 + 1 surface harmonics { YJo, Yi m , Y{* m 



m 



V 2 (l/r) = A^£l + U -^1A± = o (52) 



1,2...?}. However, in actual calculations it is often more convenient to work with a set of 

21 + 1 functions labelled by an index m running from —I to I, and to arrange things so that, 

for negative m, Yi m is also given by an expression of the form (1511) conveniently defined for 

m < 0. For this purpose, we first note the identity 

, d 2 (l/r) d 2 (l/ r) 
dr/d^ dz 2 

obtained after expressing the Laplacian operator in coordinates (£, 77, z) through the chain 
rule. The result ( 1321) suggests that we can give an operational meaning to symbols such as 
(d/drf) -1 and (d/dz)~ 2 when applied to (1/r) by means of the relations 

i d 2 = _^ 4 d_ = _ fd_\ 1 ^ 4 / d_\ 2 d_ = _ fd\ 1 

dr]di dz 2 <9£ \dr] J dz 2 \dz J <9£ \&W J 
which allow us to define V[ m for m < in terms of VJ* m | by means of (|44p . In effect, according 
to (I44p we can formally write 

/ q \ I+M / q \ -|m] / x x 

^-H-y y (-) 

and, if we substitute (d/drj)~^ by the |m|th-power of the RHS of the last operator equation 
in (133]) we obtain 



H-H = = (-4)Nl^ m| . (55) 



d\l/r) 
dz l ~\ m \di 

But it is easily seen that if we drop the absolute value sign in ( 133)) the resulting equation 



V^ m = (-l) m 4 m V^, (56) 

remains consistent with definition ( 133)) independently of the sign of m, whence (136)) can be 
freely used to relate the VJ m 's of positive and negative m. This result will be of great use 
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in Section IV on recursion relations for spherical harmonics. The corresponding relation 
for the Yim's can be now obtained from (|45p after noticing that replacing m by — m in the 
normalization constant in (lo"Tj) yields the new constant C\ _ m = C; m /4 m , or 

= C^ m r l+l V^ m = {-l) m Y? m . (57) 

The functions Yj m with —l<m<l defined by f[5"Tj) and floTl) constitute the set of standard 
spherical harmonics of degree 



IV. RECURSION RELATIONS. 



Recursion relations are of importance in calculations involving spherical harmonics, a well 
known application in Quantum Mechanics being the computation of the matrix elements of 
operators in the angular momentum representation^. In classical physics the best known 
examples are perhaps the analysis of multipolar radiation fields^^ and wave scattering 
problems solved numerically using the relatively recent Fast Multipole MethoS^-. The 
standard derivations of such recursion relations, i.e., using the associated Legendre functions 
and spherical polar coordinates, often get rather tricky and involved, yielding messy final 
expressions^. In this section we indicate an apparently new procedure to systematically 
derive a whole set of useful recursion relations using the Maxwell solid harmonics Vi m (r) 
defined in Section III [see (jSJ)]. 

We begin by noticing the equality 

Ql-m Qm ^ Ql-m Qm Ql-m-1 Qm / ^ 

l-mhoz— Hi— - > (58) 



where we have applied Leibnitz's rule to the derivatives with respect to z of (z/r). Since 
r = \/£rj + z 2 , substituting the relation 

9_ (z_\ = __iz_ = £_d_ fl\ 

dr)\r) 2r 3 / 2 2dz\r) 1 1 

into the first summand on the RHS of (1581) we obtain 

Ql—m am / 1 \ P Ql—m+1 Qm—l / \\ Ql—m— I Qm f \\ 



Q%i— m Qifj m \r J 2 Qz^~ rnJr ^ dTj m ~^ \f J dz^~ m ~^ dTj m \t 

or, in terms of the VJ m 's defined in 



2zV lm = £V Z)m _! - 2(1 - m)VU, m , (61) 
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which is our first recursion relation. A corresponding one for rjVi m can be obtained by 
replacing m by — m in (l6Tj) and taking complex conjugates using (156]) . This yields, after 
replacing m by m — 1 in the resulting equation, 

2r)Vim = -zVi^i -(l + m- (62) 

Also, adding (IBTj) multiplied by z to (|62|) multiplied by £ we find the relation 

2rV Zm = -2(7 - m)zVi- 1>m - (l + m- l)£Vj_i,m-i- (63) 

Replacing Z by Z + 1 in ( l63l) . substituting zVz m from and, finally, replacing m by m + 1 we 
obtain 

(2/ + l)£VJ m = -2r 2 V l+hm+1 + 2(l-m)(l-m- l)VJ_ 1)m+1 . (64) 

Observe that the factor £ multiplying Vi m in fl64l) depends on the angular coordinates 9 and 
if while none of the factors multiplying the functions V on the RHS does so. Therefore, 
recursion relations such as fl64l) are very convenient to compute surface integrals (projec- 
tions) of the product of the function appearing in the LHS times Maxwell's solid or surface 
harmonics, since then we can use the orthogonality properties of the latter with those on 
the RHS. Analogous expression for 7]Vi m can be obtained by replacing m by — m in ( 1641) and 
taking complex conjugates using (1561) . which yields 

2(2/ + l)r]V lm = r 2 V l+hm ^ - (I + m)(l + m - l)VU,m-i, (65) 

The corresponding expression for zVi m , is obtained by replacing m by m — 1 in fl64l) and 
substituting into flBTl) : 

(21 + l)zV lm = -rV i+1 , m -(l-m)(l + m)VU, m , (66) 

In some circumstances, it is convenient to express a derivative with respect to either 
z, 7] or £ of the product of VJ m times a function of r, say B(r), in terms of functions 
Vim multiplied by coefficients which only depend on r. This is the case, for example, in 
the computation of matrix elements of the quantum momentum operator in the angular 
momentum representation^. Since r = a/£?7 + z 2 , applying the chain rule we obtain 

|- (BV lm ) = B'zV lm /r + BV l+1>m , (67) 
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^ (BV lm ) = B'£V lm /(2r) + BV l+ i, m+1 , (68) 

^ (BV lm ) = B' V V lm /(2r) - BV l+1 , m ^/4, (69) 
where we have written B' = (dB/dr) and used the relations 

dv lm d l - m+1 d m f\\ dv lm _ d l+1 -^ d m+1 (\\ _ TA 



Using the recursion relations (161"|) - (l6"6"j) in (I6"7|) - (l6"9"j) we readily obtain the sought expressions. 

It can be checked that the recursion relations here obtained reproduce those derived in 
a different way by Kramers^ (Ch.IV p. 178). For this, it is necessary to replace in all of 
them Z by / — 1, to make B(r) = r l b(r), and to express them in terms of the (unnormalized) 
surface harmonics V\ m used by Kramers: 

9 m /1 

V lm = {-I) 1 -. ; - } r l+1 V lm . (71) 

(I — my. [I + my. 

Nevertheless, it seems to us that our derivations and expressions are simpler and more sys- 
tematic than Kramer's due to we use solid harmonics, the Vim's, instead of surface harmonics. 

Finally, of special interest in Quantum Mechanics are the ladder properties of the operator 
rx V- which is proportional to the angular momentum operator- when applied to the spherical 
harmonics Yi m . Notice that, since r x Vr - *^ 1 - 1 = 0, its effect on the Yi m 's will be the same 
as that on the functions V/ m . In order to analyze the operator r x V in a coordinate- free 
manner, it is convenient at this point to introduce the so-called spherical basis formed by 
the vector e z together with the vectors 



eg = (e z — ie y ) /y/2 and e v = (e x + ie y ) / V2 — et. (72) 
It is immediately seen that these vectors satisfy the orthogonality relations 

eg ■ eg = e,, • e,, = and ■ e v = 1 (73) 

and 

x = ie z , e v x e z = ie v , e z x eg = ieg. (74) 

Note that ( 1731) imply that eg and e,, are null vectors. In terms of the spherical vectors, the 
cartesian vectors e^ and e^ are given by 

e x = (eg + e,,) / \pl and e y = i (eg — e,^) / Vz, (75) 
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so that any vector u can be written as 

u = u x e x + u y e y + u z e z = wgeg + u v e v + u z e z , (76) 
with the components «g and u v given by 

u ■ e n = = (u x + iuy) I \Pl and u ■ eg = u v = (u x — iu y ) / \l~2. (77) 
In particular, for the position vector we have 

r = ^/V2 + T]e v /V2 + ze z , (78) 

and the gradient operator is 

^ d d d r- d r- d d 
V = e x — + e v — + e z — = V2eg — + V2e v — + e z — , (79) 
ox ay az orj at, az 



where we have used (1751) together with the relations between the partial derivatives with 
respect to the cartesian and the complex coordinates - obtained via the chain rule using ( 143]) 



d d d ^ d . d . d 
dx <9£ dr) dy <9£ dt] 

Using the recursion relations and the properties of the spherical basis derived above, it 
is easy to show that the operator e,^ • (r x V) rises the index m of functions VJ m (m = 
—I, —I + 1, — 1,/) on which it acts, while the operator eg • (r x V) lowers it. In effect, 
taking into account that 



e^ x r = ^e,, x eg/v2 + ze v x e 2 = —i£e z /v2 + ize v (81) 

and that e^ • V = y/2d/dr] [see ( T79]) ]. we find 

e^ • (r x Wi m ) = (e v x r) ■ W lm = -i£/\/2V l+hm + V2izV l+1>m+1 , (82) 

or, using recursion relation flSTj) . 

■ (r x VV lm ) = -V2i(l - m)V l>m+1 . (83) 

The corresponding property for the operator eg • (r x V) can be obtained by replacing m by 
— m in (1831) and taking the complex conjugate of the resulting expression. Thus relations 
(J56D and (EU) lead to 

eg ■ (r x VV lm ) = + m)V hm ^/V2. (84) 
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In terms of the spherical harmonics (145 D equations (|83j) and (!84j) read 

(r x VY lm ) = iy/(l-m)(l + m + l)Yi tm+1 /V2. (85) 



e r (r x Vy /m ) = V( z + m)(/ - m + l)y l)m _ a /V2. (86) 

Relations ( 185]) and (186]) are usually expressed in terms of spherical polar coordinates^, 6 
and <p. To do this, we first obtain the cartesian components of the operator 

r x V = e^— - 87 
at? sm 6/ ay? 

by recalling that = — sin ^e^. + cos <pe y and eg = e v x e r which, in turn, imply e x ■ eg = 
(e x x e v ) ■ e r = cos^cos^ and e y ■ eg = sin^cos^. Thus 

d d 
r x V = — (e x sin <p — e y cos tp) — — (e s cos cp + e y sin <p) cot , (88) 

and taking into account that e v ■ e x = 2~ 1 / 2 and • e y = i2~ 1 ^ 2 ( I85p - (I86I) can be written in 
compact form as 



,i^ f ± d_ + . cQt d_^ = ^( lTm){ i ±m + l)Yl 



mil- 



(89) 



V. LEGENDRE'S ADDITION THEOREM. 

Maxwell's integral theorem leads to a very simple derivation of the well known - and 
physically very useful - Legendre's addition theorem^. In effect, if we restrict e-^ in (138]) to 
be a real, unit vector then Hi(ei) = Y[(e\) - where Y\ = Hijr 1 - and, after dividing (138]) by 
a 1 , we obtain 

dSY.Z^ ^-^ Y^), (90) 

which must hold for any surface harmonic Y\. Therefore, if we take any set of 21 + 1 
orthonormal surface harmonics of degree I, { Xi n (r) : f dVLXi n Xi n , = 5 nn i (n = 1, ...21 + 1)} 
, the coefficients in the expansion 

21+1 

Z l (r) = J2A i X ln (r) (91) 

n=l 

are given by 

A n = j dnX* n Z t = ^fW jei), (92) 
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where we have used (!90|) with a = 1. Thus 

3(r) - 2, + i EW4(r), (93) 

n=l 

and using the relation fl35|) to express Z\ = (— l)*i!Pj(ei ■ e r ) we finally obtain the celebrated 
Legendre's addition theorem: 

4vr 21+1 

Pl(ei ■ e r ) = ^x; n ( ei )X fa (e r ), (94) 

n=l 

where we have taken into account that Xi n (r) = Xi n (e r ) since the X; n 's are surface har- 
monics. Notice that ( 1941 reproduces the familiar form of the addition theorem^ for the 
particular case in which the set Xi n (n = 0, ...21 + 1) is chosen to be the set of standard 
spherical harmonics, Yi m (r) (m = —I, .., 0, ... + I), introduced in the Section III. 



VI. SURFACE HARMONICS EXPANSIONS AND THE COMPLETENESS RE- 
LATION. 

Let us consider an arbitrarily thin sphere of radius a charged with a given surface charge 
density distribution f(9, ip). If r' = ae' r denotes the generic point on the sphere, the electric 
potential at any point r of space is 

iS- ! ^4. (95) 
r — r' 

a I I 

Inserting the expansion (J34I) into ( J95l) we obtain 

oo ; „ 

^ = E / dS ' P ^'r ■ *r)W> iff) (96) 

i=o r Jr '= a 

for r > a and, on using the expansion analogous to ( 1341) but with the roles of r and r' 
interchanged - i.e., the expansion valid for r < r' - we obtain 

oo i „ 

^) = E 4n / dS ' P ^'r ■ *r)f(e', (97) 
1=0 Jr '= a 

for r < a. Using f )96|) and f )97|) . the jump condition for the normal derivative of the potential 
at any point of the surface r = ae r (9,{p) , 4nf(9,{p) = [—d(f> / dr] a a - , writes 

00 21 4- 1 f 

w> ^ = E J dn ' Pi « ■ e ^f( e '> v0» ( 98 ) 



1=0 
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where we have set dS' = a 2 dQ'. 

Observe that if we replace each Pi(e' r ■ e r ) in the RHS of fl98|) by its expression as the 
second term of f)35p and carry out the integrals in the primed quantities, we obtain an 
expansion of f(6,ip) in terms of (Maxwell's) surface harmonics. Therefore, any finite, phys- 
ically admissible, function defined on the unit sphere admits an expansion in surface har- 
monics. The completeness of any system of orthonormal surface harmonics of the form 
[Xi n (e r ) : I = 0, 1,2... andn = 1,2. ..21 + 1] now follows by combining the results fl98|) and 
flUD as 

oo 21+1 „ 

f(B, V) = J2J2 dn 'f( e '> v')XL«)Xm(er). (99) 

l=o n=l J 

Since dQ! = sin 9'd(p'd9' , equality of both sides in ( 1991) demands the so-called completeness 
relation^: 

sin v 

1 = 71=1 

If in ( 1991) and fllOOf) we replace the generic X; n 's by the standard F; m 's we obtain the 
familiar results^ 1 ^ 

OO I p 

f{6,<p) = E E A i™Yim(e r ) -> A lm = / dn'f(9',p')YL«) (101) 

1=0 m=-l J 

and 

i=0 m=-Z 

Note that if / is independent of ip only the Aim's with m = survive in ( 11 Oil) (since 
Ym ~ e im</3 ), and therefore in this particular case we can write the expansion in terms of 
Legendre's polynomials. In effect, using ( 15T|) and replacing e' r by e z in (I35I) . we find that 
the functions Y/o P/(cos#) are proportional: 



/2TTT (-1)V +1 #(l/r) /2J+T P , m nn „, 

equation ( 1101 p can be written as 

f(0) = y)B,fl(cos0) 5, = — — / d9'sme'f{B')Pi{cose'), (104) 
where we have used the result 

[«r*»m«-<r) = saTI / = STT' (105) 
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As a first illustration of surface harmonics expansions, let us consider the electrostatic 
field due a point charge q located at a point R external to an arbitrarily thin conducting 
sphere of radius a and grounded at zero potential. This problem illustrates well the relation 
between the method of images^ and that of surface harmonics expansions 18 . The potential 
at any point r external to the sphere is 

<K r ) = r^TT + f dS'j^jr, (106) 



R 

where cr(r') is the a priori unknown surface charge distribution on the sphere. To de- 
termine it, we substitute the expansions a = J2i m ' 7 imYi rn (e' r ), l r ~ r'l -1 = 47r/(2Z + 
1) T,im r ' l /r l+1 Yim(er)Y lm (e' r ) and |r - R^ 1 = An/ (21 + 1) r l / R l+1 Y lm (e r )Y lm (e R ) into 
( I106P and carry out the surface integral using the orthogonality properties of the spherical 
harmonics, which yields: 

Anq ^ rl \ r / , v / 47r ^ a ' +2 



T J^f+T Yim( e r)Yi m (e R ) + p^i^m^m(e r ). (107) 



21 . 

Im Im 



Since = at the sphere's surface, particularizing (11071) for r = a we obtain 



aim = -^(^) l+1 Y lm (e R ), (108) 



and (11071) can be written as: 



<p(r) = I 57 , , : }. — M—Yim{e r )Yi m (e R ) = — - r — . (109) 

r — R 21 + 1 r' +i r — R r — e R a 2 /R\ 

ii i m i i i 

Equation (11091) expresses the solution to the original problem for potential at any point 
external to the sphere as the superposition of those due to q and an image (fictitious) charge 
of magnitude (a/R)q located inside the sphere at the point e R a 2 /R. Of course, if instead of 
a point charge we have a charge density distribution p(R) within a region V external to the 
sphere, we must replace q by p(R)d 3 R in dl09[) and integrate over V. 

As a second illustration, let us compute the electrostatic interaction energy of two charge 
density distributions p\ and p2 occupying non-overlapping regions V% and V2, respectively. 
Taking two points 0\ and O2 within V\ and V2, respectively, and denoting by ri (r-2) the 
position vectors with respect to 0\ (O2) of a generic point of V\ (V2), the interaction energy 
is 

E= f f dVM ^ 1 ^ 2 ] , (110) 
Jv 1 Jv 2 |r + r 2 -ri| 
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where r is the position vector of 2 with respect to 0\. Since | r 1 1 and [1*2 1 are less than |r| 
we can expand |r + r 2 — ri| _1 independently in ri and in r 2 by means of the double Taylor 

1 Yyw ^T^T \ <-) 



r + r 2 - ri fT y hi r 

lite 



But 



(-l)h( ei .V) h {l/r) PiM-er) 4tt J* d h {l/r) 



K l " *r) _. ™ \^ r Y* (r \ - K 1 ' 

h+1 2/i + l ^ hmi{ l! dzh-m^^m^ 



h\ r - L . . 

m\=— (1 

where we have used the addition theorem and expressed Yj imi (e r ) in terms of VJ imi 
using (14"5]) and (jSJ). Analogously, 

(_iy 2( e 2 . V y 2 (l/ r ) 4vr * gh(l/r) 

y 2Z 2 + 1 ^ U ' 2m2 - r ^^ e2 ^^ 2 - m2 ^m 2 ' 

m,2=—te 

whence applying the operator (e 2 • V)' 2 /^' to (I112p and commuting it with the derivatives 
in the last term we obtain, on using ()113p . 

h (eg ■ ( ei • V)M _ 4¥(-l)fa v 0*(l/r) 
Z 2 ! Zx! r (2Z 2 + l)(2fc + l) ^ * imi ' 2m2 < imi1 lJ hmal 2) dz L ^d V M ' 

(114) 

where L = h + h and M = mi + m 2 . If we now substitute (1114)1 
into (lllip taking into account that Y^Af = CLM rL+1 VLM and defining = 
[4 2 7r 2 (-l)^C /imi a 2m2 ] / [(2Z 2 + l)(2h + 1)C LM ), or 

(2/ 2 + l)(2/ 1 + l)(L + M)!(L-M)! "1 1/2 



4 2 vr 2 (-l)^ 
LM (2/ 2 + l)(2Z 1 + l) 

we finally obtain 



(2L + l)(/i + m x )!(Zi - mi)!(I 2 + m 2 )!(Z 2 - m 2 )! 



;ii5) 



E= ^rr(pi)«i^ 1 LM(e r )(p 2 )i 2m2 , (116) 

Zl/2mim2 

where we have defined the multipole moment, pi m , of a charge density distribution occupying 
a region V as 



/ dV'r n P {v')Y lm {v'). (117) 



Except for a factor of v47r, expression ( 11161) was obtained by Schwinger at al.— using a 
different reasoning. 
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VII. PARTIAL WAVE EXPANSIONS AND HOBSON'S INTEGRAL THEOREM 



Partial wave expansions are of fundamental importance in scattering problems in both 
classical 1 - and quantum physics 19 . In this section we consider the expansions in any or- 
thonormal basis of surface or solid harmonics of the function e q r - which leads to the 
partial wave expansion of a plane wave, or Rayleigh's expansion 20 -, and for the translation 
operator-!- 2 - 1 e qV , where q is a constant vector. We will see that the latter expansion imme- 
diately leads to a proof of Hobson's integral theorem^ mentioned in the introduction. Let 
us first consider the case in which q is a vector of the form a = qe q , where e q is a real unit 
vector and the scalar q may be complex. Then q • r = qr cos 9, where cos 9 = e q ■ e r , and we 
can use (11041) to expand e a r in terms of Legendre polynomials in the form 

oo 

e qr = ^/z(gr)P,(cos^). (118) 

1=0 

In order to determine the functions f\ it is first convenient to rewrite f 1 1 1 8 j) using Legendre 's 
addition theorem (1941) as 



^ = £^£^(e t )A-,„(e,.), (119) 

1=0 n=l 

where the Xi n 's are the members of any given orthonormal basis of surface harmonics of 
degree I. In terms of the regular solid harmonics H\ n = r l X\ n (I119j) writes 

oo . c. / \ 21+1 



21+1 



e qr 



If we multiply (H20p by Hi n (q)dfl q and integrate over the unit sphere using ( 1251) and taking 
into account that J dil q Hi n Hp nl = q l q v J dVtXi n X^ nl = q 2l 5u>5 nn > we obtain 

l^M^ n (r) = ^H ln {v)r l q l s l {qr), (121) 

where the function si is defined in (1241) . Hence (11191) takes the form 

oo 21+1 

e qr = A^^siiqr) ^X; n (e 9 )X, n (e r ). (122) 

1=0 n=l 

In particular, for a plane wave one has q = ik, where k is a real vector, qr = ikr and, 
since sinh(i£r) = isin(&r), it is customary to rewrite (I122p in terms of the spherical Bessel 
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functior& of order / defined as j/(C) = ( l si(i() or, using P 



m = (-D'C ( ^U^'l (123) 



Then the expansion (11221) yields the well known Rayleigh's^ expansion for a plane wave: 

oo 

eik . r = £ (2Z + i)iij^ kr ) Pl ^ k . e r ), (124) 
z=o 

where we have replaced the sum over n in fl 1 2 2 [) using Legendre's theorem ( 194)) . 
On the other hand, fl 1 2 2 [) writes in terms of the solid harmonics Hi n as 

oo 21+1 

e *' = 4nJ2 *i(<F) £ iffjq^r), (125) 

Z=0 n=l 

which is a good starting point to obtain the expansion for the translation operator- e q V . In 
effect, by formally replacing r in fl 1 2 5 [) by the operator V and qr = q^Jr ■ r by gVV 2 (see 
Section II) we obtain 

oo 2Z+1 

e «-v = 47r^J]^n(q)^n(V)5,(gVV^), (126) 

Z=0 n=l 

where we have used the fact that the operators Hi n (V) and si(q\/V^) obviously commute. 

As a first illustration of (I126p . let us consider the expansion of a function of the form 
0(|r — r'|). If, for concreteness, we assume r > r', Taylor's expansion around r = can be 
written in the well known form 0(|r — r'|) = e~ r ' v 0(r), whence 

oo 21+1 

<M|r-r'|) =47rX)E(-l) , flto('0 [^(V)^(r'Vv^)] 0(r), (127) 

i=0 n=l 

where we have taken into account that if z * n (r') = (— l) l H* n (r') An important particular case 
of (I127P is that of a spherical wave 2 -^ due to a point source located at r', e* fc ' r ~ r ''/|r — r'|, 
where k is the wave number. In this case 

( p ikr\ h.2 p ikr 
— )=- — , (128) 
r J r 

whence si(r'VV*) (e ikr /r) = si^'yT^) (e ikr /r) = si(ikr') (e ikr /r). Thus ([122]) can be 
written as 

oo 2Z+1 



ife|r— r'| 00 / ikr\ 

^- = 4^^(-l) 1 Si ( 8 k')^(r')4(V) (M (129) 



Z=0 n=l 
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or, using ( 120]) . 

p ik\r-r>\ ™ 2 J±} /I J\ l /P ikr \ 

^^-V'tWH^WHU*) (~) ( — )■ (130) 

1=0 n=l ' / V / 



It is customary to write (11301) in terms of the spherical Hankel function of the first kind and 
of order / defined as 

(i)^_(-l)'C fl d\ l fe< 



ft < «> s ^ [m) {-) • (131) 

in which case, using the definition ji(kr') = k l r' l si(ikr'), substituting H* n {Y')Hi n {y) = 
r' l r l X£ n (e' r )Xi n (e r ) and taking account that for r < r' an analysis entirely analogous to 
the one carried out above leads to an expression of the same form as fl 1 3 [) but with the 
roles of r and r' interchanged, the result (11301) can be expressed in the more familiar form 

ifc|r-r'| 00 n=2l+l 

^— = 4mkY, E Ukr^ikr^X^XU^). (132) 
l r r I 1=0 j=o 

The result (11321) constitutes the general expression for the expansion of a spherical wave in 
any basis of orthonormal surface harmonics. 

As a second application of (I126P we will provide an apparently new proof of Hobson's 
theorem^ for the integral over a sphere of radius R of the product of a surface harmonic, 
Yfc(r), and any function F(r) which is analytic for r < R. In effect, by replacing in fl 1 26[) 
vector a by r and a 2 by R 2 , Taylor's expansion around the origin for F(r) can be written 
as, 

oo 21+1 

F(r) = e r V °F = 4vr £ £ Ff n (r) [^ n (V)5,(i?Vv^)l F, (133) 

L Jo 

Z=0 n=l 

where subscript o means that the derivatives of F must be particularized at the origin. Now, 
if we normalize the surface harmonic by dividing it by \/C - where C = J dilY^Y^ -, then 
Yk/VC can always be taken as one of the elements of the generic basis Xi n , in which case 
r k Yk/\J~C replaces one of the Hi n in ( 11330 . The rest of the basis elements could be obtained, 
for example, via Gram-Schmidt orthogonalization. Therefore, since Yj. is orthogonal to all 
solid harmonics in H133[) except to H k /yC and J dSY k Hl = R k+2 C, we obtain 

/ dSY k F(r) = 47rR k+2 \s k (RV^)H k (V)] F, (134) 

Jr=R 1 J ° 

which is Hobson's integral theorem. 
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VIII. THE ROTATION MATRIX FOR THE SPHERICAL HARMONICS. 



It is well known 2 ^* 2 ^ that if a given rotation brings the position vector r into the vector 
r' - the image of r-, then the value at r' of any spherical harmonic of degree I, Y$ m '(r'), is a 
linear combination of the values at r of all 21 + 1 spherical harmonics of the set of degree 
/, i.e. Yj m /(r') = Y^ m =-i ^rn'mXimi^) ( m ' — where the elements of the rotation 

matrix, D l m , m , are uniquely determined by the rotation. In this section we will carry out 
an analysis of the rotational properties of the spherical harmonics based on the invariance 
under rotations of the scalar product, that is, if r' and b' are the images of vectors r and b, 
then (b • r) = (b' • r') or, since r = r', (b • e r )' = (b' • e' r ) . Thus if b is a null vector, so is 
its image b' (b 2 = b' 2 = 0) and (b • r)' and (b' • r/ are solid harmonics whose values at e r . 
and e' r can be expanded in terms of the Y/ m (e r ) and Yi m (e' r ) ( (m = —I, ...I), respectively. 
Then the elements of the rotation matrix are obtained after equating corresponding terms 
in both expansions. 

According to ( 11011) we can write 

oo I „ 

( b ■ = E E / dQ " ( b • O l VmiO, ( 135 ) 

1=0 m=-l 

where the double prime denotes the generic point on the unit sphere. Since, on the unit 
sphere (b ■ r)' = (b ■ e r ) 1 , the integrals on the RHS of (I135[) can be computed using Maxwell's 
theorem (13T1) with 



*.(b.r)' and Y , sYit= J»±l (-D' + ^^ a-(l/r) 

V 4tt ^(l + m)\(l-m)\dz l - m di m 

which yields 



/ 



«>»<! „»\>v*<„»\ I in (_i)-. 2 - 9 (b ■ r) 

da (b • e„) Y lm (e r ) = V( , + m)!(t _ m) , ^-„ as „ ("7) 

In order to compute the RHS of fl 1 3 Tj) . we express b in the spherical basis, 

b = 6 5 e 5 + b v e v + b z e z , (138) 

and dot it with r = (£e^ + i]e v )/\/2 + ze z , which yields 

r i l 1 /'2~- J / 2 

(b-r)'= (b^ + b v C)/V2 + zb z =E /,'_,iui fa + ^ H i'; ] . (139) 
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Thus, since the only term in the binomial expansion of (b • r) which contributes to the 
derivative in the RHS of f!137j) is that containing the product 2:' -m £ m - whose coefficient is 
readily found from (I139p to be l\2~ m l 2 b™b l ~ m /[(I — m)!m!]- we have 

d l (b • r) ; 
where we have used the fact that 



^| 2~m/2 yntf-m _ ^ r^j2-m, jl-m^(l-m)/2^(l + m)/2 (140) 



b 2 = -> b z = 2b^b rj or b • e z = 2b ■ e v b ■ e ? . (141) 
Substituting (I140p into (11371) we finally obtain the expansion 



/ i m=l - m , (/-m)/2,(«+m)/2 v / v. 

(b ■ r/rf = W iu[j^ T LjL h Ylm(er \ (142) 
VaZ + l^ v/(/ + m)!(/-m)! 1 J 

which was introduced by Kramers^ and, later, by Schwinger— as a generating function for 
defining the spherical harmonics on the LHS. Here, we have preferred to proceed instead by 
expanding (b • r/r) 1 into the already known [see (1511) ] Yi m J s and, then, finding the expansion 
coefficients using Maxwell's integral theorem. 

The invariance under rotations of (I142j) makes it a suitable expression to analyze the rota- 
tional transformation properties of the spherical harmonics. In effect, let there be given two 
orthonormal, righthanded sets of vectors, (e x ,e y ,e z = e x x e y ) and (e' x ,e' y ,e' z = e' x x e^J, 
and consider the rotation about the origin of coordinates defined as that bringing the vec- 
tor e^. into e' x and the vector e y into e' y . Then, any vector which rotates rigidly with the 
unprimed system will be brought by the rotation into another vector - the image - that has 
the same components with respect to the primed system as has the former with respect to 
the primed one. Hence if with respect to the spherical basis associated with the unprimed 
set a vector b has components b%, b v and b z and its image, b', has components b'g, U q and b' z 
, we have 

b' = b'^ + b' v e v + b' z e z = & f eJ + h n e' v + b z e' z , (143) 

where e^, e' v and e' z are the vectors of the spherical basis associated with the primed set. On 
taking the scalar product of both sides of the last equality in (I143j) by e' and e^, respectively, 
we obtain 

k = • e ? + b 'v e 'v ■ e v + & H • e *> ( 144 ) 
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and 

b v = &Je£ ■ e ? + b' v e'^ ■ e v + b' z e'^ ■ e z . (145) 
Now, if b is a null vector, so is its image under rotation, b', whence 



K = iyjWp v (146) 
and, since and are null vectors, (I141j) yields 



e' v - e z = i^J 2e' v - e v e' v ■ e c and • e z = -iyj 2e^ • e„ • e c , (147) 

where we have taken different branches for the square root in order to satisfy the condition 
• e 2 = (e^ • e z ) . Introducing ( I146j) and f!147j) into (11441) and f!145j) it is easily seen that 
the transformation laws for the £ and i] components of a null vector can be written as the 
perfect squares 

and 

h n = (b'l /2 v/^V^ + K 1/2 V^V^) 2 • (149) 

Thus, defining the quantities s = b 1 ^ 2 , t = b^ 2 , c = ■sj e' v ■ and d = e'^ ■ , we can 
write (I148I) - (I149I) more compactly as 

s = cs' - d*t', (150) 

t = ds' + c*t'. (151) 
It is easy to show that (I150I) - (I151I) is an unitary transformation with unit determinant - 

1/2 1/2 

i.e. it belongs to the group SU2-, so that the quantities b^ and b v transform as the 
components of a two- spinoi^^-. In effect, from the definitions of c and d, and according to 



(11471) we have cdc*d* = \Jg'^ ■ e v e'^ ■ e^e'^ ■ e^e'^ ■ e n = e' v ■ g z g'^ ■ e z /2 and, therefore, 

(cc* + dd*) 2 = e' v ■ e'^e'^ ■ + e' f • e^e' v ■ e v + ■ ■ e z — e' n ■ I ■ e'^ — 1, (152) 

where 

I G X G X — GyGy -\- G Z G Z G^G^ — G^G^ -\- G Z G Z (153) 

is the unit dyadic. Since we assume that no reflections are involved in the transformation of 
the coordinate system we take the positive root in HI 53[) and obtain cc* + dd* = 1, whence 
the transformation fl 1 5 [) - f 1 1 5 1 [) is unimodular. 
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Now, if the images under the rotation of the null vector b and of the position vector r 
are b' and r', respectively, an expansion analogous to (j!42p must hold for (b' • r'/r') 1 where 
&£, b v and e r are replaced by b'g, b' v and e' r . Thus, equating both expressions we obtain 



' i m 's' l - m 't' l+m 'Y lm/ (e' r ) i m s l ~ m t l+m Yi m (e r 



To find Yi m i(e' r ) in terms of the Y5 m (e r )'s we just have to substitute the transformation 
fll50l) - fll5ip into the RHS of (I154p and, for each m, find the term containing s ,l - m 'f l + m . 
The factor multiplying s' l ~ m ' t' l+m ' in the m-th summand can be computed as 

J Ql—m' Ql+m' 

{l- m ')\{l + m ')\ds' l - m ' dt' l + m ' ^ C$ ' ~ ^ m(rfS ' + cY ) /+m ] ( S ',t')=(0,0) ' ( 155 ) 
or, applying Leibnitz rule to the (Z + m')-th partial derivative with repect to t' of the product 
in (I155P and evaluating the result at t' = 0, 

1 d l ~ rn ' (I -m)l(l + m )\(-iYd* k c* l+rn '- k c l - m - k d m ~ m ' +k s a - m ' 
(I - m')\ ds n - m ' f^ Q k\(l + m' -k)\(l-m- k)\(m - m! + k)\ ' ^ 156 ^ 

with the convention that the factorial of a negative number is infinity. Carrying out now the 
derivative with respect to s' in (11561) and substituting the result in (I154p we readily obtain 

i 

Y lm ,{9'^') = D^Je'^e^Y^e,^), (157) 

m=—l 

where (9' , ip') and (9, (p) are the spherical coordinates of the points e' r and e' r , respectively, 
and we have set 



l+m' 



k=0 



1 I - w! i 


II 


'l + m'\ 


\ m — m! + k j 




I * J 



(-l}k d *k c *l+m>-k c l-m-k d m-n,>--!, 




(158) 



with 

c = = ^eVe*, and d = • e € = ^/e'* • e*. (159) 

Expression fl 158ft yields the elements of the rotation matrix, D l m , m (e' v , e v ), for the rota- 
tion of the coordinate system defined as that bringing the spherical vector into e' v - or, 
equivalently, the vector e x into e' x and the vector e y into e' y -. Notice that (1 1 5 8 [) - fl 1 5 9[) are in- 
dependent of the particular parametrization chosen to describe the rotation, and that (I159P 
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provides with a direct geometrical meaning for parameters c and d of the SU2 transforma- 
tion (I150p - (I15ip . In particular, it is easy to show that for the usual Euler angles employed 
in Quantum Mechanics 22 to parametrize a rotation, namely: a) a rotation by an angle a 
around the z-axis, b) a rotation by an angle /3 around the new y-axis and c) a rotation by 
an angle 7 around the new z-axis, one obtains (by expressing the vectors [e^, e'] in terms of 
[e x , e y ] and the Euler angles) 

(e'„ • e;) 1 / 2 = e - 4(Q+7)/2 cos(/3/2), (e'J • e* v ) 1 ' 2 = e l ^~ a)/2 sin(/3/2), (160) 

and (I158p goes into the standard form of the rotation matrix elements^. 

Finally, note that we have considered here only the case of rotation matrices of integer 
rank. However, in the theory of quantum angular momentum^ 2 , both integer and half- 
integer ranks are, of course, of importance. A very complete presentation of compact, 
parametrization-free representations of finite rotation matrices of arbitrary rank using in- 
variant tensor forms and spinor operators can be found in the works by Manakov et al.— 7 — , 
where they also give some applications such as photo-emission by polarized atoms or reduc- 
tion formulae for some tensor constructions characterizing photoprocesses. 



IX. GAUNT INTEGRALS AND WIGNER COEFFICIENTS. 

In this section we will apply Hobson's integral theorem (I134p to obtain a generating func- 
tion for integrals of products of three spherical harmonics - the so-called Gaunt integrals 2 ^- 
and provide their expressions in terms of the well known Wigner coefficients^ 2 -. For this pur- 
pose, let us first use ( 1142 j) to express the regular solid harmonic associated to the spherical 
harmonic Y lm , H lm = r l Y lm , as 



Hi m (r) 



(21 + + m)l(l - m)\ 



(b ■ rV 



, i+ ' (161) 



47r(-iy+ m 2 2 '/! 2 

where the subscript on the RHS means the coefficient of s l ~ m t l+m = ^ _m )/ 2 5^ +m )/ 2 j n the 
binomial expansion of (b • r) 1 . Thus we have 



(2/ fe + 1)(2Z C + l)(2l d + l)fx b \ii c \fx d \u b ^cW- 



r UJ ■ r2 hm b - r2 l c ma- r2 ldm d ~ ]l 47r( — 1) A/ 2 2L /! 2 /! 2 /! 2 1 PbPcPd> 



where 



I = J dQ (b • r) lb (c • r)' c (d ■ r) ld , 



(163) 
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and, to simplify notation, we have written 

fi = l-m, u = l + m, p = s l - m t l+m , 

L = l b + l c + l d and M = m h + m c + m d . (164) 

In order to compute (I163p we apply Hobson's integral theorem f)134|) with R = 1, Yk — 
Hk = (b • r) lb and F(r) = (c • r) ic (d • r)' d , whence 



47T 



s Z6 (VV2)(b- V)' 6 [(c-r) /c (d-r) id ] . (165) 



Since (c ■ r) /c (d • r)' d is a homogeneous polynomial of degree l c + Id, it is clear that the RHS 
of (11651) vanishes unless A& = Z c + Id — h is zero or a positive even integer. Therefore, the 
only term in the power series (1241 for Si b (V 2 ) which contributes to (11651) is that containing 
(V 2 ) Afc / 2 , which corresponds to j = L/2, where L = l c + l d + h- Thus, taking into account 
d2U), (HS5]) yields 

1 = ^ V2 ) Ab/2 ( b ■ V )"]o K c ■ r )"( d ■ r)' d ] • ( 166 ) 

Now, since V 2 (c ■ r)' c = V 2 (d ■ r)' d = we have 

V 2 [(c ■ r) lc {d ■ r) ld ] = 2V(c • r)'° • V(d ■ r) ld = 2l c l d c- d (c • r) h ~ 1 (d ■ r) 1 ^ 1 , (167) 

and repeating the process A&/2 times we easily obtain 



(V 2 )^/ 2 [(c • r)Hd ■ r) 1 *} = ^|]|^^ (c • r) A ^ 2 (d • r)^ 2 , 

where we have written X c = l d + h — l c and A^ = l b + l c — h- Thus, 

/ = ^ L/2 (L/2)U c U d l r „ , Ad/2 Ac/2] 

(L + l)!(V2)!(A c /2)!(A,/2)! Llb Vj J « LlC rj ld T) J ' 



(168) 



(169) 



Also, Leibnitz's rule yields 

h 

I Ac/21 = 

./!(/, ./)! 



(b ■ V)' 6 [(c ■ r) Ad / 2 (d ■ r) Ac / 2 ] = ^ fc " [(b • V) j (c • r) Ad / 2 ] (b • V)^'(d • r) Ac / 2 , 



3=0 

(170) 

and the only term in the series in the LHS that gives non zero contribution at r = corre- 
sponds to j = Xd/2, which also makes h — j = A c /2. Thus 

[(b • V) ib ] o [(c • r) Ad / 2 (d • r) Ac / 2 ] = l b \(b ■ c) Ad / 2 (b • d) Ac / 2 , (171) 
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and inserting (j!7ip and (I169P into (11661) we finally obtain 

An2 L / 2 (L/2)\l b \l c U d l (c ■ d) A ^ 2 (b ■ c) A ^/ 2 (b • d) A ^ 2 
~ (L + l)!(A b /2)!(A c /2)!(A d /2)! ' 

Notice that (I172p vanishes if is odd or if any of the parameters At, A c or A^ 
the values of l c and ^,7^0 only if the value If, is one of the set 

h + ld, lc + ld-2, ■■■ \lc-l d \- (173) 

The reader familiar we group theory might have noticed that this set of integers is the same 
as that in the Clebsch-Gordan series^ for the decomposition into irreducible representations 
of the tensor product of two irreducible representations. Also, observe that the foregoing 
procedure based on Hobson's integral theorem can be extended, in principle, to compute 
integrals of products of any number of simple solid harmonics (b • r)' 6 (c • r)' c ...(q • r) z « but, 
of course, the result of applying the operator [(V 2 )' ic+ " i<r ~' 6 " 2 (b ■ V)' 6 ] o gets more involved 
as the number of factors increases. 

In order to compute I PbPcPd we must find the coefficient of the term s^ b t u b b s^ c t v c c s^ d t v d d in the 
expansion of the product (c ■ d) Ai, / 2 (b ■ c) Ad / 2 (b ■ d) Ac / 2 where X b = l c +ld—h, K = ld+h~L 

1 /2 1 /2 

and X b = h + l c — Id- Note that if we write b = b^e^ + b v e v + b z e z , Sb = b^ , i& = b v and 
b z = i\/2s&t& (b null vector), with corresponding definitions for s c , t c , s d and t d , we have 

b ■ c = s\tl + t\s\ - 2s b t b s c t c = (s b t c - t b s c ) 2 (174) 

and, analogously, 

b ■ d = (s b t d - t b s d ) 2 and c ■ d = (s c t d - t c s d ) 2 , (175) 

whence 

(c • d) A ^ 2 (b • c) Ad / 2 (b • d) A ^ 2 = (s c t d - t c s d ) Xb (s b t c - t b s c ) Xd (s b t d - t b s d ) x \ (176) 
By binomially expanding (s^ c — tbS c ) Xd , 

( Sb t c - t bSc ) x « = ( ~ 1)Ad "^.f 4 t x r j tisc d - j , (177) 

(Ad-jJIj! 

we see that the only term in the binomial expansion of (sbt d — tbS d ) Xc which multiplied by 
s{t b d ~^ in (11771) yields s^t^ = pb is that containing s b t b c ~ k with k = — j + fi b which, in 



(172) 
Thus, given 
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addition, also makes A c — k + Xd — j = X c + Xd — //& = Ub- Therefore, 

*6*«J - Mi ^ s^c - * 6 s c Ad ^ = ]T \ 1 d c c d d . 178 

pb b b {*d- jjhKto - j)K u b - Xd + jy- 

Also, the term containing s^t Xb ~ k t^s d q '~ k with k = j — X d + fi c is the only one in the 
binomial expansion of (s c t d — t c Sd) Xb which can yield s£ c t£ e s^ d t^ d = p c pd when multi- 
plied by t J c s x ^H^ b ~ J s ~ Xd+Ub+: > in (11781) . It is clear that such a value of k also makes 
Xb — k+j = Xb + Xd — p c = v c and, since we must also have that v d = k + p b —j = p b — Xd + p c , 
the m-indexes must satisfy the condition 

M = m b + m c + m d = 0. (179) 

That relation (1 179ft must be satisfied for (II 62ft to be non vanishing becomes obvious if we 
express the integral on the LHS in spherical coordinates, since it assures the independence of 
the azimuthal coordinate of the integrand. Observe that (I179P also implies that the exponent 
of Sd is Xb — k — Xd + v h + j = Xb + v h — ji c = fi d , as it should be. Therefore, 



J PbPcpd 



[(c-d^/^b-c^OD-d)^/ 2 ] 

(-iy^ +m ^x b \x c \x d \ (i8o) 

^J\(X d -j)\^ b -j)\(iy b -X d + j)\(j + ^ c -X d )\(iy c -jy: 

where the sum in (j!80p is over all integers j for which the factorials all have nonnegative 
arguments. Inserting (11801) into the RHS of (11721) and the later into fl 1 6 2 [) we finally obtain, 
after dividing by R L+2 , 



1 h Ic Id ' 


II 


' h Ic Id J 


1 I 




\m h m c m d J 



cMY lbmb Y lcmc Yi dmd = \I K \ n I II I (181) 

where the integral in the LHS vanishes unless 

m b + m c + m d = 0, l h + l c > Id, l c + k> h and U + h> k, (182) 
and we introduced the standard Wigner 3j-symbols 29 defined as 



k Ic Id = y (-l) lb+l ^ v /A{lbJ c Jd)Vp b W^dW'ycW. (lg3) 



where 



Hh,Ud)^^0^, (184) 
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and 




(-l)^/yA(M c ,/,)(L/2)! 
(A 6 /2)!(A c /2)!(A d /2)! 



(185) 



X. CONCLUSIONS. 

We have presented the coordinate-free approach to spherical harmonics based on the 
contributions by Maxwell, Thomson and Tait. It has been shown that many results of the 
theory of spherical harmonics needed in physical applications can straightforwardly and sys- 
tematically derived from the simple properties of elementary solid harmonics (b • r)' ( b null 
vector), Maxwell's harmonics and Maxwell's integral theorem. In particular, we have pro- 
vided simple and apparently new proofs of known results such as Maxwell's integral theorem 
and Legendre's addition theorem, the latter being obtained in terms of a general basis of 
orthonormal surface harmonics. Also, it has been shown how the derivation of recursions 
relations can be made simpler and more systematic using Maxwell's solid V/-harmonics than 
- as is conventionally done - using surface harmonics and spherical coordinates. Surface 
harmonics expansion have been discussed and their relation with the method of images il- 
lustrated and, also, they have been used to efficiently compute the interaction energy of 
non-overlapping charge distributions. We have reviewed in a unified manner partial wave 
expansions and provided a seemingly new proof of Hobson's integral theorem. A procedure 
based on compact vector methods has been given to express the elements of the rotation 
matrix in a form which is independent of the parametrization of the rotation. Finally, we 
have shown how Hobson's integral theorem can be used to find integrals involving products 
of three (or more) spherical harmonics. 

Appendix A: Derivation of the expression for Yi m in spherical coordinates. 

The Vrharmonic VJ m defined by (jSJ) can be easily expressed in spherical coordinates 
through the familiar Legendre polynomials as follows. First, consider the well known Leg- 
endre's expansion 




(Al) 
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where Pj(fi) is the Legendre polynomial of degree j, 

P( AZ ) = J_^L( /i 2_ 1 v (A2) 
Taking m times the derivative with respect to \x in both sides of flAlj) we get 

I 1 v A m P 1 x P 

* 1 y^gj-w g j _ 1 V s i+m fA3l 

(1 - 2su + s 2 ) m+1 /2 (2m + 1)!!^ rfu m (2m + 1)!!^ d/i m ' 1 ; 

where we have taken into account that d m Pj/dn m = for j < m. Making s = h/r and 
// = z/r = e z • r/r, with r = (x 2 + y 2 + z 2 ) 1 ^ 2 , we can write ( 1A3I) as 

1 r" 2m A ^ (i m P. 



[x 2 + y 2 + {z - / i )2]™+i/2 ( 2m + i)!! ^ r i+i ^ 

Now, it is readily seen that we can compute g l - m (i/ r 2m + ls )/Q z l - m as (_iy- m times the 
(/ — m)-th derivative of (1A4j) with respect to /i evaluated at h = 0, or 

= f iV-m ^~ m ) ! r -G+™+i)gjg (A5) 
9z^r w (2m + 1)!! (i/x m ' 

Therefore, 



' 2 m r' +m+1 dri m 2 m r l+1 drt m 

Substituting (1A6j) into (!5Tl) and using (IA2I) we obtain the familiar-^ expression for Y\ 



which is valid also for m < 0. 
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